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Abstract. We study the efficiency of numerically simulated galaxy clusters for producing strong gravitational lensing events in 
different dark-energy cosmologies with constant and time- variable equation of state and we compare it with "standard" ACDM 
and OCDM models. Our main results are that: (1) the expected abundance of gravitational arcs with large length-to- width ratio 
depends on the equation of state of dark energy at the epoch of formation of the halo; (2) as found in previous studies, the strong 
lensing cross section of galaxy clusters is very sensitive to dynamical processes like mergers, whose impact however is different 
for different cosmologies, being stronger for models in which halos are less concentrated. As expected, the largest differences 
in the lensing optical depth occur at moderate and high redshift. 



1. Introduction 

Evidence is mounting that the Universe is spatially flat, has 
low matter density and is dominated by some form of dark 
energy, acting as a repulsive gravitational force and responsi- 
ble for the present phase of accelerated cosmic expansion, (e.g . 
[ Riess et alJl 19981 IPerimutter et alJI 19991 lHawkins et alJl2003l 
Ispereel et alJ 120031 iTeemark et alJ l2004 iRiess et alJ l2004 . 
The dark energy generalises the concept of cosmological con- 
stant, admitting dynamics and fluctuations. Such a generali- 
sation is necessary attempting to reconcile its value with the 
known energy scales of particle physics. The minimal exten- 
sion of a cosmological constan t is a self-in teracting scalar 
field (Quintessence, see lPeebles & Ratrall2002l and references 
therein), al ready considered before the evidence fo r cosmic ac- 
celeration ( Wetterich 1988; Ratra & Peeble sll988l) . 

Structure grows earlier in dark-energy than in 
models, which makes dark- 



matter haloes more concentrated 
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Following earlier analytic work ( 


Bartelmann et al. 


2003) 



we study here with cluster models numerically simulated in 
eight different cosmologies how earlier structure growth and 
denser halo cores affect strong lensing by galaxy clusters. 
Our motivation is three-fold. First, despite recent attempts 
(Wambsg anss et alJl2004t iDalal et alJ l2003). it remains to be 
clarified whether the strong-lensing efficiency of real clusters 
is reproduced by numerical cluster models or not. Earlier 



studies found that the expected strong-lensing efficiency of 
clusters in cosmological-constant models appears to fall short 
by an ord er of magnitude of repr oducing the observed number 
of arcs (Bartelma nn et alJ Il998t) . We thus wish to quantify 
the effect on strong lensing of the increased halo core density 
in dark-energy models. Second, earlier structure growth may 
have interesting implications for strong cluster lensing in 
view of the recently-detected clusters at redshifts near unity in 
which strongl y-lensed arcs were found (Hasinger et al. 1998; 
Gladders et al. 2003; Zaritskv & Gonzalez 2003; lLidman et all 
2004j). Th ird, we found in a recent study with high time 
resolution ( Torri et ai1l2004 that the efficiency of a cluster for 
strong lensing can substantially be increased during merger 
processes because of the enhanced tidal field and surface-mass 
density during different phases of the merger process. Shifting 
structure growth to earlier redshifts also allows major merg- 
ers at higher redshifts and thus potentially much increased 
arc-formation probabilities. 



With these motivations in mind, w e use the 17 c l uster- sized 
haloes produced for and described in IDolag et alJ §004) and 
investigate their strong-lensing properties. Section 2 describes 
the eight cosmological models used, Sect. 3 the cluster sample, 
and Sect. 4 the lensing simulations. Results are presented in 
Sect. 5 and summarised in Sect. 6. 
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2. Cosmological models 

We compare in this paper the efficiency for producing strong 
lensing events of the sampl e of numerically s imulated clus- 
ters previously described in Dol ag et all J2004I) . We give here 
a short description of both the cosmological models and the 
methods used for obtaining the cluster sample and refer the 
reader to that paper for further details. 

All clusters in the sample were simulated in different cos- 
mologies with constant and time-variable equation of state. 
These are an open Cold Dark Matter (OCDM) and four 
spatially flat cosmological models, namely a cosmological- 
constant (ACDM) model, a dark energy model with constant 
equation of state (DECDM), and two quintessence models 
with markedly different dynamical properties, one with inverse 
power-law Ratra-Peebles potential (RP, see Peebles & Ratra 
2002, and ref erences t herein) and one with SUGRA potential 
(SUGRA, see iBrax & MartirJ l2000. and references therein). 
For the latter three we have used two different normalisations 
of the perturbation power spectrum, to be discussed later. 

In all cases, the matter density parameter today is Qq = 
0.3. In the flat cosmologies, the remaining 70% of the critical 
density is assigned to the dark energy at present (Oa = 0.7. 
The remaining cosmological parameters are h = 0.7, Qb h 2 
0.022. Gaussian density fluctuations with scale-invariant power 
spectrum, and no gravitational waves. 

The normalisation of the perturbation power spectrum, 
which we set conventionally through the rms density fluctu- 
ation level erg within spheres of 8/z _1 Mpc radius, has been 
chosen as erg = 0.9 for all the cosmological models in or- 
der to match the observed abundance of galaxy clusters. Since 
the power spectrum normalisation is most important here, af- 
fecting the structure statistics and therefore the lensing power, 
we also normalise the DECDM, RP and SUGRA models to 
the observed level of an isotropics in the C osmic Microwave 
Background (CMB) fe.g. lBennett et alfeooll) . In this case, the 
erg is generally slightly smaller because of the Integrated Sachs 
Wolfe (ISW) effect affecting the large scale CMB anisotropics: 
the reason is the larger time interval in which the cosmic equa- 
tion of state changes in dark energy scenarios with respect to 
the cosmological constant case. This enhances the dynamics 
of the gravitational potentials and raises the ISW power, vis- 
ible on t he large angle tail of C MB anisotropy power spec- 
trum ("seelBartelman n"et alJl2002l and references therein). We 
take 0.86, 0.82 and 0.76 for DECDM, RP and SUGRA, re- 
spectively. Note that the reduction is quite important, since 
we selected our dark energy models in order to highlight the 
differences among each other and the cosmological constant. 
If the equation of state is smaller than the cases we consider 
here, the ISW effect and the reduction to the power spec- 
trum normalisation are smaller. These numbers as well as all 
the inputs from the linear evolution of cosmological perturba- 
tions to the N-body procedure described later are computed 
using our dark energy oriented cosm ological code, based on 
CMBfast (Seliak & Zaldarriaga 1996). The code, originally in- 
troduced bv lPerrotta & Baccigalupildl999l) . describes the dark 
energy as a scalar field, numerically integrating its background 
and linear perturbation dynamics for the most relevant track- 



ing quintessence scenarios, capable to converge to the present 
phase of cosmic acceleration fro m a wide ra nge of initial condi- 
tions in the early universe (Baccigalup i et alJ2000l) . This is per- 
formed in the general framework of scalar-tensor cosmologies, 
in which the d ark energy may be directly linked with gravity 
(tPerrotta et al.l200ol) . The outputs of this machinery have been 
used for other work, allowing in particular to constrain the dark 
energy equ ation of state with the m odern CMB anisotropy ob- 
servations (Baccigalu pi et alj20 02). 

The dark energy cosmologies are characterised by their 
equation of state parameter w(z), describing the ratio between 
pressure and density of the Dark Energy. For cosmic expansion 
to be accelerated today, the w parameter must be < -1/3. In 
the limit of constant w, the dark energy density evolves with 
redshift proportionally to (1 + z) 3(I+ll) . That is the case for both 
the ACDM model, whose w = — 1 yields a constant dark en- 
ergy density, and the DECDM model, which we choose to have 
w = -0.6. It is useful to note that from the point of view of the 
Friedmann equation, the curvature term in the OCDM model 
beha ves as a dark energy co mponent with constant w = -1/3 
(see lBartelmann et all2002t and references therein). 

For the RP and SUGRA models, w is time-dependent. 
The dark energy is consistently described by means of the 
quintessence scalar field </>. The RP and SUGRA quintessence 
potentials are given by 



M 4+a e 
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■4+ffn 



SUGRA - 



■exp(47rG<^ 2 ) , 



(1) 



<j) a Q (j) a Q 

respectively . The exponential in V sur.RA reflects super-gravity 
corrections jBrax & Martinl2 000) and induces large variations 
in the equation of state with respect to the RP case at the epoch 
of structure formation, even if w is the same at present in the 
two models. Ipolag et all (|2004) found indeed that the cluster 
concentration reflects the dark energy behaviour at the time of 
collapse, a promising feature in order to reduce the degener- 
acy between different dark energy models yielding the same 
amount of acceleration today. In this work we check how this 
feature affects the statistics of strongly lensed arcs in clusters. 
Specifically, the tracking RP and SUGRA quintessence mod- 
els considered here have wq = -0.83 today, with czq = -0.6 
and ffg = -6.7, respectively. In the tracking regime, this yields 
wsugra — -0.23 and wrp -0.77. 

3. Cluster sample 

The cluster models used in this study were simulated us- 
ing a new version of the cosmological code GADGET 

j * 1 II i 

(Springel et al. 2001), specifically adapted to simulate the for- 
mation and evolution of cos mic structures in t he dark energy 
models described earlier (see lDolag et all2004l) . 

Each cluster in our sample was obtained using the so- 
called re-simulation technique, which consists in re-simulating 
at higher resolution a patch of a pre-exist ing large-scale cos- 
mological simulation ( Torme n et alJ 119971) . For this work, we 
selected as a parent simulati on an N-body run w i th 512 3 parti- 
cles in a box of479/T 1 kpc dYoshida et all200lUjenkins et alJ 
2001 ). Its background cosmological model is spatially flat with 
£2 m ,o = 0.3 and Q A = 0.7 at the final epoch, identified with 
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redshift zero. The particle mass was 6.8 x 10 10 /z 1 M Q , and 
the gravitational softening was chosen as 30h~ l kpc. From the 
output of this simulation at z = 0, we selected ten spherical 
regions of radius between 5 and 10 Mpc, each containing 
either one or a pair of dark matter haloes with mass exceeding 
3 x 10 14 M Q . We thus obtain a sample of 17 cluster-sized 
objects. 

Each of these regions was re-sampled to build new initial 
conditions with on average 10 6 dark-matter particles. The ini- 
tial conditions from the ACDM cosmology were adapted to 
all the dark energy cosmologies studied here as described in 
Dola g et alJ fe004 . The mass resolution of the re-simulation 
ranges from 2 x 10 9 to 6 x 10 9 h~ x M Q per dark-matter particle 
so as to have each cluster consisting of approximately the same 
number of particles. The gravitational softening is reduced to 
a 5 h' 1 kpc cubic spline smoothing for all particles in the high- 
resolution region. 

We found that our 17 clusters contain on average Ny ~ 
200, 000 dark matter particles within their virial radii. The cor- 
responding virial masses range between My = 3.1 X 10 14 to 
1.7 x 10 15 h~ l M e . 

4. Lensing simulations 

Ray-tracing simulations are then carried out with each of the 17 
clusters in all eight cosmological models. The technique used 
in th is study was described in detail in several ear lier papers 
(e.g. lBartelmann et alJl995lMeneghetti et alJ200ol) . 

For all clusters in each cosmological model, we produce 
^Vsnap = 52 snapshots at different redshift between and 1. 
These are equidistant in time, leading to a time resolution of 
~ 100 Myr, which allows the effects of dynamical processes on 
the cluster c ross sections for s trong lensing to be resolved, as 
discussed in lTorri et alJ yJXM). 

For each snapshot, we select a cube of 3 /; _1 Mpc comoving 
side length, centred on the halo centre and containing the high- 
density region of the cluster. The particles in this cube are used 
for producing a three-dimensional density field, by interpolat- 
ing their position on a grid of 256 3 cells using the Triangular 
Shaped Cloud method ( Hockney & Eastwood 1988). Then, we 
project the three-dimensional density field along the coordinate 
axes, obtaining three surface density maps £y, used as lens 
planes in the following lensing simulations. The total number 
of lensing simulations is therefore 17x52x3x8 = 21216, re- 
quiring approximately the same number of computation hours 
on an IBM-SP4 cluster located at the computer centre of the 
Max Planck Society in Garching. 

The lensing simulations are performed by tracing a bundle 
of 2048 x 2048 light rays through a regular grid, covering the 
central quarter of the lens plane. This choice is driven by the 
necessity of studying in detail the central region of the clusters, 
where critical curves form, taking into account the contribution 
from the surrounding mass distribution to the deflection angle 
of each ray. 

Deflection angle s on the ray grid are c omputed using the 
method described in Meneg hetti et We first define a 

grid of 128 x 128 "test" rays, for each of which the deflection 



angle is calculated by directly summing the contributions from 
all cells on the surface density map Sy, 

4G v x Kk - Xjj 
^ = J ]x^x~J ' (2) 

where A is the area of one pixel on the surface density map 
and Xh,k and X/j are the positions on the lens plane of the "test" 
ray (h, k) and of the sur face density element (z, j). Following 
IWambsganss et all dl998l) . we avoid the divergence when the 
distance between a light ray and the density grid-point is zero 
by shifting the "test" ray grid by half-cells in both direc- 
tions with respect to the grid on which the surface density is 
given. We then determine the deflection angle of each of the 
2048 x 2048 light rays by bi-cubic interpolation between the 
four nearest test rays. 

The position y of each ray on the source plane is calculated 
by applying the lens equation. If y and x are the angular po- 
sitions of source and image from an arbitrarily defined optical 
axis passing through the observer and perpendicular to the lens 
and source planes, this is written as 

y - x - -^a(x) , (3) 

where D\ s and D s are the angular diameter distances between 
the lens and the source plane, and between the observer and the 
source plane, respectively. 

Then, a large number of sources is distributed on the source 
plane. We place this plane at redshift z$ = 1. Keeping all 
sources at the same redshift is an approximation justified for 
the purposes of the present case study, but t he recent d etections 
of arcs in high-reds hift clusters dZaritskv & Gonzalez! 120031 
Idladders et afll2003l) indicate that more detailed simulations 
will have to account for a wide source redshift distribution. 

The sources are elliptical with axis ratios randomly drawn 
from [0.5, 1]. Their equivalent diameter (the diameter of the 
circle enclosing the same area of the source) is r e = 1". These 
are distributed on a region on the source plane corresponding 
to one quarter of the field of view where rays are traced. As 
in our earlier studies, we adopt an adaptive refinement tech- 
nique when placing sources on their plane. We first start with 
a coarse distribution of 32 x 32 sources and then increase the 
source number density towards the high-magnification regions 
of the source plane by adding sources on sub-grids whose reso- 
lution is increased towards the lens caustics. This increases the 
probability of producing long arcs and thus the numerical effi- 
ciency of the method. In order to compensate for this artificial 
source-density enhancement, we assign a statistical weight to 
each image for the following statistical analysis which is pro- 
portional to the area of the sub-grid cell on which the source 
was placed. 

By collecting rays whose positions on the source plane 
lay within any single source, we reconstruct the images of 
background galaxies and measure their length and width. 
Our technique for imag e detection and c lassification was de- 
scrib ed in detail by Bartelmann & Weiss ( 1994j) and a dopted 
bv iBartelmann et alJ dl998l). iMeneehetti et alJ feOOOl l200ll 
l2003albh and Torri et al. (2004). It results in a catalogue of sim- 
ulated images which is subsequently analysed statistically. 
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5. Results 

In this Section, we discuss how the lensing properties of the 
clusters in our sample differ in the different cosmological mod- 
els. 



5.1. Critical curves and caustics 

We start by considering the critical curves, along which 
strongly magnified and highly distorted images form. These 
are the lines on the lens plane where the determinant of the 
Jacobian matrix of the lens mapping, 



Ahk{X) = = Ohk - -7, 



(4) 



is zero. Since the local magnification is the inverse of the 
Jacobian determinant, the critical points are ideally charac- 
terised by infinite magnification. 

Through the lens equation, critical lines are mapped onto 
the caustics on the source plane. Sources lying close to the 
caustics therefore have images strongly magnified in the radial 
or tangential directions, i.e. gravitational arcs. The probability 
of one cluster to produce highly distorted images is then re- 
lated to the size and the shape of the critical curves and of the 
caustics. 

In Fig^we show some examples of how the critical curves 
(and of the corresponding caustics) of one of the clusters in 
our sample change among different cosmological models. The 
halo has a virial mass ranging between 4.3 x 10 14 M Q and 
5.2 x 10 14 M G at redshift z ~ 0.3. Considering models char- 
acterised by the same normalisation of the power spectrum 
of the primordial density perturbations, it is clearly seen that 
the two models with dynamical dark energy, i.e. the RP and 
SUGRA models, and the DECDM model with constant equa- 
tion of state w = -0.6, interpolate between the ACDM and the 
OCDM model. In particular, the critical curves in the SUGRA 
model have sizes comparable to those in the DECDM model. 

However, it is important to note that the size of the criti- 
cal curves and of the caustics in dark energy models may be 
smaller when the CMB normalisation is used, because of the 
enhanced ISW effect with respect to the ACDM case. The 
larger the required reduction of erg, is, the more the critical 
curves shrink; in the case considered here, the critical curves 
tend to have sizes comparable to those in the ACDM model or 
smaller. 

5.2. Lensing cross sections 

In earlier work, the length-to-width ratio of gravitational arcs 
was found to be very sensitive to many intrinsic properties of 
the lenses which depend on cosmology. In fact, numerical sim- 
ulations showed that order-of-magnitude differences are ex- 
pected in the number of arcs on the whole sky with length- 
to-width ratio exceedin g a given t hreshold in differe nt cosmo- 
logical mod els (Bartelmann & Weisslll994l iBartelmann et alJ 
1 19951 Il998l) . In particular, ten times more arcs are expected 
in an OCDM than in a ACDM model. Modelling clusters using 
the NFW density profile, which depends on cosmology in con- 



trast to the singular isothermal profile, such a high sensitivity 
of the abundance of long and thin arcs on cosmology has been 
confirmed even using an ana lytic approach (Meneg hetti et alJ 
b003blEartelmann et all2003l) . 

Therefore, aiming at evaluating the differences between the 
strong-lensing efficiency of clusters in different cosmological 
models with dark energy, we focus on the statistical distribu- 
tions of the arc length-to- width ratios. 

The efficiency of one galaxy cluster for producing arcs 
with a given property can be quantified by means of its lens- 
ing cross section. By definition, this is the area on the source 
plane where a source must be placed in order to be imaged 
as an arc with that property. As explained in the previous sec- 
tions, we artificially increase the density of background sources 
close to the cluster caustics by adopting adaptive grid refine- 
ment. Each source is taken to represent a fraction of the source 
plane. We assign to each source and all of its images a statis- 
tical weight w which is inversely proportional to the squared 
resolution of the sub-grid on which the corresponding source 
was placed. The cells of the sub-grid with the highest resolu- 
tion have area A, and the sources placed on its grid points are 
given a statistical weight of unity. The absolute lensing cross 
sections are then determined by counting the statistical weights 
of the sources whose images have a length-to-width ratio ex- 
ceeding a threshold (L/W) m ; n . If a source has multiple images 
with (L/W) > (L/W) m j n , we multiply its statistical weight by 
the number of such images. Therefore, the lensing cross section 
is 



„„ = A ^ WiWirii 



(5) 



where W, is unity if the i-th source has images with (L/W) > 
(L/W) m i n and zero otherwise, n,- is the number of images of 
the i-th source satisfying the required condition, and w, is the 
statistical weight of the source. 

5.2.1 . Averaged cross sections 

We now discuss the averaged cross sections of our cluster sam- 
ple. In Fig. [2] we show the lensing cross sections for long and 
thin arcs, averaged over all 17 clusters, as function of the lens 
redshift for all eight cosmological models. For each cluster, we 
use its three projections along the coordinate axes for measur- 
ing the cross section. Therefore, the curves displayed in Fig. [3] 
result from averaging over 5 1 different curves in each cosmo- 
logical model. 

As expected, the lensing cross sections reflect the differ- 
ences in the concentration of dark matter halos in different cos- 
mological models. This has b een fully discussed in our previ- 
ous paper JPolag etalJ E004). It is shown there that the clus- 
ter concentration depends on the dark energy equation of state 
at the cluster formation redshift z co n through the linear growth 
factor D + (z co n). Assuming the same normalisation of the power 
spectrum, the lensing cross sections for the OCDM and the 
ACDM models differ by roughly a factor of four, independent 
of the minimal length-to-width ratio of the arcs, and the cross 
sections for the other cosmological models with dark energy 
interpolate between them. Despite the equation of state of the 
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Fig. 1. Critical curves (upper panels) and caustics (bottom panels) of one of the clusters in our sample in all eight cosmological 
models. The side length of each panel is 50" 




Fig. 2. Averaged lensing cross sections for arcs with length-to-width ratio larger than 5 (left panel), 7.5 (central panel) and 
10 (right panel) of our cluster sample as function of the lens redshift. Sources are kept at redshift z s — 1. Different panels 
refer to different cosmological models. Solid curves show the results for cosmologies with erg, = 0.9. Dashed curves show the 
correspondent results when cr% is reduced for taking into account the increasing Sachs-Wolfe effect affecting the large scales in 
theCMB. 



dark energy today is the same for the RP and SUGRA models, 
their lensing cross sections differ substantially at higher red- 
shift. 

For models with the CMB normalisation of the power spec- 
trum, we find lensing cross sections smaller by even more than 
one order of magnitude compared to the OCDM model. In fact, 
when the normalisation is reduced because of the Integrated 
Sachs Wolfe (ISW) effect affecting the large-scale CMB 
aniso tropics in the cosmologies we consider (B artelmann et alJ 
not only the formation epoch of our simulated clusters is 
delayed, but also their evolution up to redshift zero is changed. 
For example, clusters in the RP and in the SUGRA models have 
masses at redshift zero which are smaller by roughly 20% and 
30% respectively compared to the ACDM model. 

The cluster sample is still too small for the averaged cross 
section to be a smooth function of redshift. In fact, the curves 
exhibit strong peaks which are caused by merger events aris- 
ing in single clusters. It has recently been shown that during 
such events, on timescales of some hundreds of Myr, the cluster 
effici ency for strong lensing is strongly enhanced dTorri et alJ 
120041) . due to the combined effects of the increasing shear and 
convergence. We will discuss the properties of these peaks in 



the next subsection, when the lensing cross sections of one in- 
dividual cluster will be shown. 

5.2.2. Individual clusters 

In Fig. |2 several peaks can be found in at least all the curves 
obtained in cosmological models with the same normalisation 
of the power spectrum. Indeed, in this case, we are comparing 
the same clusters in different cosmological models and there- 
fore the same dynamical processes cause the cluster lensing 
efficiencies to peak. However the position and the amplitude 
of the peaks is strongly dependent on the cosmological model. 
This becomes clearer by looking at the lensing cross sections 
of one individual cluster. 

In the left panel of Fig|3]we show the lensing cross sec- 
tion for arcs with length-to-width ratio larger than 7.5 of one of 
the most massive clusters in our sample (~ 10 15 h~ x M ). Solid 
curves refer to models with <x 8 = 0.9. Since the cluster forms 
at different epochs in different cosmological models, the same 
merger events occur at different redshifts, and consequently the 
peaks are shifted in time. Initial conditions of our clusters are 
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Fig. 3. Left panel: lensing cross sections for arcs with length-to-width ratio larger than 7.5 of one of the most massive clusters 
in our sample as a function of the lens redshift in different cosmological models. As in Fig. |3 sources are kept at redshift 
Zs = 1. Solid curves show the results for cosmologies with erg = 0.9. Dashed curves show the corresponding results for <xg 
reduced according to CMB normalisation. Central and right panels: redshift evolution of the derivative of the mass with respect 
to redshift and of the velocity dispersion of the same cluster as in the left panel. Solid and dashed lines refer to the different 
normalisations of the power spectrum. 



built such that single halos reach the same evolutionary state at 
redshift zero in all cosmological models. Therefore, the delay 
between dynamical processes in different cosmologies shrinks 
as the redshift decreases. 

The offset between the merger histories is clear from the 
central and the right panel of Fig. [3] In the central panel, the 
redshift evolution of the redshift derivative of the halo mass 
within the virial region, defined as the region enclosing 200 
times the mean background density, is shown. The mass grows 
towards redshift zero as the cluster accretes material from its 
surroundings. Sudden peaks in the curves indicate that large 
clumps o f matter are enter ing into the virial region. As dis- 
cussed in iTorri et all J2004I) . the typical time-scale for an in- 
falling substructure to reach the peri-centre of its orbit is then 
~ 1 Gyr, when the passage of the substructure produces a peak 
in the particle velocity dispersion, as shown in the right panel 
of Fig. [3] All the most pronounced peaks in the curves repre- 
senting the lensing cross sections correspond to some peak in 
the velocity dispersion, which demonstrates our previous asser- 
tion about the relation between merger events and maxima in 
the cluster lensing efficiency. 

By comparing the left and the right panels of Fig.[3J we note 
that we cannot find for all the peaks in the velocity dispersion 
a corresponding peak in the lensing cross sections. The reason 
is that peaks at high and low redshift are suppressed by the ge- 
ometrical lensing efficiency. Indeed, the lensing cross sections 
drop when the lenses are too close to the observer or to the 
sources. 

It is remarkable that even if the lens is intrinsically not very 
efficient for lensing, the passage of massive substructures close 
to the cluster centre can strongly enhance its ability to produce 
gravitational arcs. For example, in the ACDM model the lens- 
ing cross section is larger by roughly one order of magnitude 
while a merger is occurring. Similar increments in the ampli- 
tude of the lensing cross section have been found bv lTorri et alJ 



l2004> . This suggests that lensing events are likely to be tran- 
sient events, characterising peculiar epochs during cluster evo- 
lution. 

When the CMB normalisation of the power spectrum is 
adopted, the comparison between the same clusters in different 
cosmological models is less straightforward. Indeed, changing 
erg strongly affects the merging history of the cluster, as shown 
by the dotted lines in Fig. [3] 

We also note that in both the RP and the SUGRA models 
with CMB normalisation, even when the cluster appears to be 
dynamically perturbed, the lens still does not become particu- 
larly efficient for strong lensing. We verified that in these cases, 
the substructures producing the peaks in the velocity dispersion 
are orbiting at large distance from the cluster centre. Indeed, the 
merger events producing the peaks in the lensing cross sections 
in the cluster-normalised models are here delayed in time due 
to the reduced <x 8 . 

5.3. Optical depth 

We can now estimate the probability for arc formation of a pop- 
ulation of clusters. Given a mass function, dn(M, z)/dz, assum- 
ing that the typical lensing cross sections for lenses of mass M 
are well represented by the lensing cross sections of our sam- 
ple <t(M, z), the optical depth r(z s ) is defined as the sum of the 
lensing cross sections of each lens between the observer and 
the sources, divided by the area of the source plane. Following 
this definition, 



T(Zs) 



1 



AnDi 



dz(l +zY 



dV 



dz 



f°° dn 

dM 

Jo dM 



<t(M,z), (6) 



where z s is the source redshift, V is the cosmic volume, and 
the factor (1 + z) 3 accounts for the fact that the cluster number 
density is defined per comoving volume. The former equation 
applies under the assumption that lensing cross sections do not 
overlap. 
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Fig. 4. Left panel: differential optical depth for arcs with L/W > 7.5 and sources at z s = 1 as a function of the lens redshift. 
Solid curves show the results for cosmologies with cr 8 = 0.9. Dashed curves show the corresponding results when <x 8 is reduced 
according to CMB normalisation. Right panel: optical depth for arcs with L/W > 7.5 and sources at z s = 1 as a function of the 
normalisation of the power spectrum. 



Since we do not know cr(M, z) for a continuous range of 
masses, but only for a discrete cluster sample M,, we calculate 
the optical depth as follows: 



T(Zs) = 



1 



4ttD- jo 

n-l 

X 



dz(l +zf 



dV 



dz 



Y J \[(r(M i ,z) + cr(M i+u z)] 

i 

JM; 



(7) 



^ Mi+1 dn(M,z) 
x I ... dM 

J Mi 



dM 



where n is the total number of clusters in our sample. This cor- 
responds to attaching to each halo of mass M in the interval 
[Mj,M i+ i] a lensing cross section cr(M,z) = l/2[cr(M,,z) + 

(T{M M ,Z)1 

The integrand in Eq. 7 is shown in the left panel of Fig.|4] 
for different cosmologies. The differential optical depth for arcs 
with length-to-width ratio larger than 7.5 and for sources at 
redshift z s = 1 is larger in those cosmological models where 
lenses form earlier and are thus more concentrated. Moreover, 
the contribution to the total optical depth comes from clusters 
in a wider redshift range in these cosmologies. For example, in 
the high redshift tail, the curves drop to zero at z ~ 0.65 and 
z ~ 0.8 in the ACDM model and in the OCDM models, re- 
spectively. Moreover, at z ~ 0.6, the differential optical depth 
is still close to its maximum in the OCDM model, while it is 
decreased below 30% in the ACDM. Other cosmologies, like 
the RP, the DECDM and the SUGRA with cluster abundance 
normalisation, interpolate between these models, while we ob- 



tain substantially smaller optical depths by adopting the CMB 
normalisation of the power spectrum. 

The total optical depth as a function of the normalisation of 
the power spectrum is shown in the right panel of Fig. |4] Again, 
we find that when cluster normalisation is used, dark energy 
models interpolate between ACDM and OCDM. When the nor- 
malisation is reduced according to the amplitude of large-scale 
temperature fluctuations in the CMB, the optical depth also de- 
creases as shown by the solid lines in the figure. 

The total optical depth changes by a factor of two between 
the ACDM and the OCDM models. This is in disag reement 
with the previous results of iBartelmann et aD Jl998l) . where 
one order of magnitude larger optical depths were found in the 
OCDM compared to the ACDM. However, we use here mod- 
ified method for estimating the optical depth. Moreover, both 
the mass and in particular the time resolutions of our simu- 
lations are considerably larger than in the earlier paper. This 
allows us to better account for the contribution to the optical 
depth from all the merger events occurring in the lenses, all 
of which enhance the cluster efficiency for producing arcs, as 
shown in the previous sections. The impact of mergers appears 
to be different in different cosmological models, being larger 
in those cosmologies where lenses are less concentrated. As a 
consequence of that, the difference in the number of arcs which 
is expected to be observed in a ACDM and in a OCDM cos- 
mological models is reduced when mergers are accounted for. 
Indeed, resampling our cross sections for re producing the same 
time resolution of Bartelma nn et alJ {l998), we obtain an opti- 
cal depth larger by a factor of eight in the ACDM compared t o 
the OCDM model, in agreement with (Bartelma nn et al.jT9 981. 
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d [h- 1 Mpc] 

Fig. 5. Cross sections for arcs with L/W > 7.5 as a function 
of the distance between two merging halos, modelled using a 
pseudo-elliptical NFW lens model with ellipticity e = 0.3. The 
first halo has mass M ma ; n = 7.5 x 10 14 h~ l M Q , the second M su b = 
3 x 10 14 /z _1 M o . The lens and the source planes are placed at 
Z\ = 0.3 and z s = 2, respectively. Different line types refer to 
different cosmological models. 

This interpretation is also supported by the following sim- 
ple experiment. We simulate the passage of a clump of mat- 
ter of mass M su b = 3 x 1Q 14 h~ l M Q through the centre of an 
halo of mass M main = 7.5 x IO'VMq. Both the main clus- 
ter clump and the merging substructure are modelled using the 
pseudo-ellipti cal Na v arro-Fr enk- White lens model discussed 
by iMeneghett i et alJ (|2003b), assuming an ellipticity of the 
lensing-potential contours of e — 0.3. The lens and source red- 
shifts are kept fixed at z\ = 0.3 and z s = 2, respectively. Having 
placed the main halo at the centre of the lens plane, we produce 
deflection-angle maps for several distances between the merg- 
ing clumps by moving the substructure along the x-axis. Then, 
we apply the methods discussed in the previous sections for 
determining how the lensing cross sections evolve as functions 
of the distance between the colliding halos. Results are shown 
in Fig. [5] Different line types refer to different cosmological 
models and all curves are normalised to the cross section of an 
unperturbed halo of mass = 7.5 x \0 l4 h~ l M Q . 

We note that in this analytic experiment we reproduce well 
the three peaks in the lensing cross sections which were pre- 
viousl y found in numerical lensing simulations by iTorri et alJ 
(2004). As discussed there, the first and third peaks are due to 
the increasing shear in the region between the two merging ha- 
los, while the substructure approaches the main cluster clump. 
The second peak is due to the larger convergence of the lens 
when the two clumps are perfectly aligned. As anticipated ear- 
lier, the growth of the lensing cross sections relative to the un- 
perturbed case is larger by more than 50% in the ACDM than 
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Fig. 6. Optical depth as a function of the minimal lens redshift 
contributing to the integral in Eq. 7. All cross sections of lenses 
below redshift z were neglected. 



in the OCDM model. Other cosmological models with static 
and dynamical dark energy interpolate between ACDM and 
OCDM, confirming our assertion that mergers have a stronger 
impact in cosmological models where halos are less concen- 
trated. 

Bart elmann et alJ (|2003) already estimated using analytic 
models the change of the optical depth relative to the ACDM 
model in different cosmological models with constant equation 
of state of the dark energy. The only cosmological model we 
investiga ted here, which can be directly compared to the re- 
sults of Barte lmann et alJ (|2003) is the DECDM model with 
CMB normalisation. In this particular case, we find an optical 
depth which is larger by roughly 25% compared to the ACDM 
case, while the analytic models predicted a larger increment, 
i.e. between 50% and 90%. This inconsistency is again due to 
the smaller impact of mergers in the DECDM compared to the 
ACDM model. 

As discussed earlier, we see larger differences in the dif- 
ferential optical depth among different cosmological models at 
high redshift. This is better shown in Fig. [6] where we plot the 
integrated optical depth as a function of redshift, when the con- 
tribution from lenses below redshift z is neglected. The curves 
are normalised to the ACDM model. The curves show that the 
deviations from the ACDM model become significantly larger 
at redshifts higher than ~ 0.6. Therefore, searching for arcs in 
high-redshift clusters can be very promising for constraining 
both the cluster evolution and the equation of state of dark en- 
ergy. 
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6. Summary 

We have investigated the efficiency of numerically simulated 
galaxy clusters for producing strong lensing events in different 
cosmological models with constant and time-variable equation 
of state of the dark energy, and we have compared it with "stan- 
dard" ACDM and OCDM models. Both the local abundance of 
galaxy clusters and the observed amplitude of the temperature 
fluctuations in the CMB have been used for normalising the 
power spectrum of the primordial density perturbations gener- 
ating our sample of numerically simulated clusters. In the latter 
case, the normalisation is reduced for accounting for the in- 
creasing ISW effect affecting in different ways the large-scales 
CMB anisotropics in the cosmologies we considered. When the 
local cluster-abundance normalisation is used, the initial con- 
ditions in the N-body simulations were set up such as to obtain 
halos having approximately the same evolution at redshift zero. 
Our main results can be summarised as follows: 

- the lensing cross sections for long and thin arcs averaged 
over all clusters in our sample reflect the differences in the 
concentration of dark matter halos in different cosmologi- 
cal models. Halos are more efficient for producing strong 
gravitational lensing events in those cosmological models 
where they are more concentrated. Assuming the cluster 
abundance normalisation of the power spectrum, we find 
lensing cross sections which are smaller by roughly a fac- 
tor of four in the ACDM compared to the OCDM model, 
and the other dark energy models interpolate between them. 
Moreover, despite their present equation of state is the same 
for the RP and the SUGRA models, their lensing cross 
sections differ substantially. For those dark energy mod- 
els where CMB normalisation is adopted, we find lensing 
cross sections smaller by more than one order of magnitude 
compared to the OCDM model, indicating that the power 
spectrum normalisation is a crucial aspect for assessing the 
relative behaviour of the strong-lensing effect in different 
dark energy cosmologies. 

- We find that the local maxima in the lensing efficiency are 
caused by merger events occurring in the lenses. In partic- 
ular, we verify that the peaks in the redshift evolution of 
lensing cross sections of individual clusters are shifted in 
time due to the offset between the merger histories in dif- 
ferent cosmologies. 

- The optical depth for lensing is larger in those cosmolog- 
ical models where lenses form earlier and are thus more 
concentrated. Moreover, the contribution to the total opti- 
cal depth comes from clusters in a wider redshift range in 
these cosmologies. 

- Compared to the ACDM model, the optical depth in the 
OCDM model is larger by only a factor of two, in disagree- 
ment with the results of previous other work, which found 
order-of-magnitude differences between these two models. 
However, the higher time resolution of our simulations al- 
lows us to account for the merger events occurring in the 
lenses, which we verified to have a stronger impact in those 
cosmological models where lenses are less concentrated. 

- When the same normalisation of the power spectrum is 
adopted, we find that clusters in the dark-energy models 



we investigated have optical depths within the boundaries 
set by the ACDM and the OCDM models. However, when 
the CMB normalisation is used, the optical depths become 
significantly smaller. The differences between the cosmo- 
logical models are relatively small, but they become larger 
when the contribution to the lensing optical depths from 
high redshift clusters is considered. 

Based on these results, we conclude that arc statistics is a 
potentially very powerful tool for constraining the equation of 
state of dark energy, and for investigating the dynamical evo- 
lution of galaxy clusters. In particular, strong lensing events 
at high redshift represent an important source of information 
which might be used for discriminating between dark energy 
models. Further investigations are needed now for fully under- 
standing the impact of mergers on the lensing cross sections 
and its dependence on cosmology. 
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